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Levy flights and fractional Brownian motion (fBm) have become exemplars of the heavy tailed 
jumps and long-ranged memory widely seen in physics. Natural time series frequently combine both 
effects, and linear fractional stable motion (Ifsm) is a model process of this type, combining a-stable 
jumps with a memory kernel. In contrast complex physical spatiotemporal diffusion processes where 
both the above effects compete have for many years been modelled using the fully fractional (FF) 
kinetic equation for the continuous time random walk (CTRW), with power laws in the pdfs of both 
jump size and waiting time. We derive the analogous kinetic equation for Ifsm and show that it has 
a diffusion coefficient with a power law in time rather than having a fractional time derivative like 
the CTRW. We discuss some preliminary results on the scaling of burst "sizes" and "durations" in 
Ifsm time series, with applications to modelling existing observations in space physics and elsewhere. 

PACS numbers: 89.75.k,89.75.Da 
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I. INTRODUCTION 

Fractional kinetics is finding increasingly wide applica- 
tion to physics, chemistry, biology and interdisciplinary 
complexity science [l|, H, S 0, B Q- One reason for 
this is the link between "strange" kinetics and observed 
non-Brownian anomalous diffusion, motivating the use of 
fractional dynamical models of transport processes, in- 
cluding those based on fractional calculus. Scalas et al. 
[5|] give numerous applications; we will just note space 
plasmas 01, magnetically confined laboratory plasmas 
[3) B IC J llli ) fluid turbulence ^ and the travels of dollar 
bills [12. 

Equally widespread in its application, and evolving in 
parallel with the theory of anomalous diffusion, is the the- 
ory of anomalous time series. The corresponding mod- 
els, particularly in the mathematics and statistics liter- 
ature, have often been based on stable self-similar pro- 
cesses [13, HJ, [i3- Stability here means the property 
whereby the shape of a probability density function (pdf ) 
remains unchanged under convolution to within a rescal- 
ing (c.f. chapter 4 of Mantegna & Stanley [la]) It is 
an attractive feature in modelling, particularly when one 
anticipates that a signal represents a sum of random pro- 
cesses. In particular stable self-similar processes, a devel- 
opment in the wider field of stochastic processes |17l KLSJ. 
can model two effects which are are often seen in real 
data records. The first-Mandelbrot's "Noah" effect [T^- 
describes non-Gaussian "heavy-tailed" pdfs, while the 
second-his "Joseph" effect _20] -manifests itself as long- 
ranged temporal memory. The many applications have 
included hydrology [21| , finance [l0| , magnetospheric ac- 
tivity as measured by the auroral indices 2^ 2^1 , 
solar wind quantities [23|, and solar flares 
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in-situ 



The existence of two rich, parallel, but intersecting lit- 



eratures means that it is not yet completely known which 
techniques from one will apply to a given problem in the 
other. It is for example not always clear a priori what 
type of kinetic equations will apply in a given context. 
The right class of kinetic equation for reversible micro- 
physical transport need not also be the right one for an 
evolving time series taken from a macroscopic variable. 
The problem of model choice is an important and timely 
one, both in physics and more general complexity re- 
search. Because different models can predict subtly dif- 
ferent observable scaling behaviours, distinguishing them 
may require measuring several exponents, as any individ- 
ual exponent may be identical across several models, a 
point recently emphasised by Lutz [25|. 

It is now increasingly recognised that much natural 
data is of the type that Brockmann et al. [l^ dubbed 
"ambivalent" . In such systems heavy tailed jumps and 
long-ranged temporal memory compete to determine 
whether transport is effectively superdiffusive or subd- 
iffusive. The ambivalent process they used for illustra- 
tion, and fitted to data, was the well-studied fully frac- 
tional continuous time random walk j3| which incorpo- 
rates both effects via fractional orders of the spatial and 
temporal derivatives in its kinetic equation. Zaslavsky 
et al. |26l . |27| also advocated use of the same process in 
space physics for modelling the auroral index time series. 
They explicitly contested [271] the applicability of a time 
series model [23] based on a self-similar stable process- 
linear fractional stable motion (Ifsm) [ij, [ij]-in this role. 
We note that, rather than being purely a mathemati- 
cal abstraction, Ifsm has been linked to physics via the 
propagation of activity fronts in extremal models [28|. 
A comparison of these two approaches, leading to a bet- 
ter understanding of their structure and their similarities 
and differences, thus seems to us to be highly topical. It 



will be the first of two main topics of this Paper. Al- 
though we are fully aware that the kinetic equation we 
obtain on its own cannot fully specify a non-Markovian 
process, and, importantly, will not be unique to Ifsm, we 
nonetheless believe that our comparison of the kinetic 
equations for the two paradigms is of value, particularly 
as a source of physical insight (see also sections 1 and 2 
of 0). 

To make the comparison we first briefly recap (Table 
1) the main kinetic equations corresponding to the mod- 
elling of time series by stable processes, and of anomalous 
diffusion by the continuous time random walk (CT RW) , 
respectively. In particular we highlight (following |25| l 
the difference between fractional Brownian motion (fBm) 
and the fractional time process (ftp) which has sometimes 
led to confusion, at least in the physics and complexity 
literature (e.g. [2^). We illustrate the potential value of 
this comparison with reference to a surprising gap in the 
physics literature, the absence of a kinetic equation cor- 
responding to Ifsrn, analogous to the one given for fBm 
by Wang & Lung [30] ■ We give a simple derivation by di- 
rect differentiation using the characteristic function. The 
kinetic equation can be obtained by methods as diverse 
as a transformation t"^ of the time variable in the space 
fractional diffusion equation, and a path integral [3ll |. 

Our second main topic is the potential relevance of 
Ifsm to physics as a toy model for "calibrating" diagnos- 
tics of intermittency (c.f. [34|)- As a frequent attribute 
of nonequilibrium and nonlinear systems, intermittency 
has been a particular stimulus to physicists and time se- 
ries modellers |33j . In particular the paradigm of self- 
organised criticality (SOC) [33] has been one framework 
for this, embodying the hypothesis of avalanches of activ- 
ity in nonequilibrium complex systems. We investigate 
the scaling of intermittent bursts in Ifsm, using the burst 
size and duration measures which have very often been 
used as direct diagnostics of SOC. Such measures have 
been previously studied on, among many others, mag- 
netospheric and solar wind time series [3J|. We follow 
several earlier conjectures [22, [SJ, [3^, [3y| and make sim- 
ple scaling arguments building on a result of Kearney 
& Majumdar [37| which suggest that Ifsm could indeed 
be one candidate model for observed power laws in such 
bursts. We test our arguments with numerics using the 
algorithm of Stoev and Taqqu [3^| . We confirm the ear- 
lier numerical results of Carbone et al. [35| and the more 
recent work of Rypdal and Rypdal [3^. These papers 
considered just the a — 2, fBm case, using a running 
average threshold and a fixed threshold respectively. 

However we find numerically that for Ifsm our simple 
scaling argument, while giving a good approximation to 
the dependence of the burst size exponent on the self sim- 
ilarity exponent iJ as a is progressively reduced from 2, 
becomes much less accurate as a reaches 1. We conclude 
by offering some suggestions as to the reasons for this, 
and describing future work. 
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H 


Stable process 
H ^1/a + d 


CTRW 
d! ^aH 


a=2 


H = 1/2 


BWBm 
V^p = dtp 


< a < 2 


H ^ 1/a 


oLm 
V'P = dtp 


Q = 2 


0<H<1 


fBm 
2Ht^''-^V^p = dtp 


ftp 


< Q < 2 


0<H<1 


Ifsm 


ap 

v"p = dr"p 



TABLE I: Kinetic equations for the main classes of process 
used to study anomalous time series and transport beyond the 
Bachelier Wiener Brownian paradigm. Ordinary Levy motion 
(oLm) parameterised by a stability exponent a relaxes the fi- 
nite variance assumption of the central limit theorem. The 
fractional time process (ftp), and the ambivalent process (ap) 
i.e. the fully fractional continuous time random walk, add 
a fractional derivative of order aH to the kinetic equations 
for WBm and oLm respectively. A different way of introduc- 
ing temporal memory effects is via stable selfsimilar processes 
with memory kernels, fractional Brownian motion (£Bm)and 
linear fractional stable motion (Ifsm) respectively. Although 
not fully specified by them, the stable processes nonetheless 
have kinetic equations with time dependent diffusion coeffi- 
cients. To our knowledge the (boldfaced) kinetic equation for 
Ifsm had not been derived before our preprint [2^. It was 
subsequently arrived at by path integral methods in [3l|]. In 
the self-similar stable processes the self similarity parameter 
H depends on a memory exponent d and on the stability ex- 
ponent a via H — 1/a + d. In the CTRW case by contrast 
the standard memory exponent, d' — aH. In all cases there 
is a coefficient D with appropriate dimensionality on the left 
hand side which we have set to 1, in the BWBm case this is 
simply the familiar diffusion coefficient. 



II. THE KINETIC EQUATION OF LFSM 

A. Limit theorems and stochastic processes 

In Table 1 we collect the kinetic equations for the pdf 
p = p{x, t) of some processes which have been proposed in 
the various literatures on time series analysis and anoma- 
lous diffusion. For clarity we concentrate on the simplest 
examples from the family of stable processes and from 
the CTRW. In the table the (statistical) self similarity 
exponent H is defined using dilation in time where At 
goes to AAi: 



a;(AAi) = X"x{At) 



(1) 



and the equality is in distributions. 

The fourth row of Table 1 corresponds to Bachelier- 
Wiener Brownian motion (BWBm) and the fifth row to 
the familiar diffusion equation where we have abbrevi- 
ated d/dt to dt- BWBm is of course a manifestation of 
the central limit theorem (CLT) [13, EI- The solution 
p{x, t) is of Gaussian form with width spreading as t^ " 



and its characteristic function is also a Gaussian in k 
(~ exp(— |fcpt) with stabihty exponent a = 2; the pro- 
cess has power spectrum S{f) ^ /"^. 



B. Anomalous diffusion and the extended Central 
Limit Theorem 

Similarly, the fifth row corresponds to relaxing the as- 
sumption of finite variance, by allowing a stability expo- 
nent < a < 2. The corresponding probability density 
functions (pdf) Paix, t) are the a-stable class, with power 
law tails decaying as x~^°'^^\ Following Mandelbrot we 
refer to these as "Levy flights" or ordinary Levy motion 
(oLm). 

The corresponding kinetic equation 



dpa{x,t) _ d°'pa{x,t) 



at 



d\x\°' 



(2) 



has a symmetric Reisz fractional derivative in space, 
9"/9|a;|", which in the Table is given as three di- 
mensional and abbreviated to V". The Reisz deriva- 
tive is a pseudodifferential operator with symbol — |fc|" 
and Pa{x,t) has characteristic function pa(k,t) — 
exp(— |fc|"i). Unlike the cases we now go on to dis- 
cuss, the kinetic equation for oLm is still unambiguously 
Markovian and an expression of the extended CLT. Due 
to infinite divisibility, in this specific case the pdf alone 
Pa{x, t) is enough to uniquely characterise the stochastic 
process, which we will call Za(t). 



C. Relaxing independence through temporal 

memory: fractional BroAvnian motion vs. the 

fractional time process 

The seventh row of Table 1 describes the case when 
the iid assumption is relaxed, rather than the finite vari- 
ance one. This case is more subtle than the previous 
two. Relaxing independence is one way to break the iid 
assumption and is the situation we consider. It can be 
done in several ways, we will discuss just two. 

One of the ways which has been employed in the 
CTRW formalism is to take a power law pdf of waiting 
times p{t) ^ t^'^^+"-^^[43]. This became known as the 
fractional time process (ftp, see also [25|). The waiting 
times themselves are still iid, but their infinite mean is 
assumed to be a consequence of dependence due to mi- 
croscale physics. The kinetic equation that corresponds 
to the ftp |25l . I43I |44| can be seen in the fourth column of 
the eighth row in Table 1 . We may define a temporal ex- 
ponent by d — aH. The fractional derivative in time, of 
order aH — d , corresponds physically to the power law 
in waiting times. The prime indicates that this exponent 
is not identical to the memory parameter d in the case 
of fBm or FARIMA ^^. d' runs from to 1 and is, for 
example, the same as the temporal exponent defined by 
Brockmann et al (their "a"; our a is their "/3"). In all 



the following cases D is no longer the Brownian diffusion 
constant but simply ensures dimensional correctness in 
a given equation. Note that we do not include the term 
describing the power law decay of the initial value here 
or in subsequent CTRW equations (it is retained and dis- 
cussed in [43], see their eqn. 40). 

Another way to relax independence is to introduce 
global long range dependence, as pioneered by Mandel- 
brot & Van Ness [21|. They used a self-affine process 
with a memory kernel, originally due to Kolmogorov and 
called by them "fractional Brownian motion" . Contradic- 
tory statements exist in the physics literature concerning 
the equivalent kinetic equation for fBm corresponding to 
that for ftp. It has somtimcs been asserted [1, 23] that 
the equations are the same, while ftp has sometimes been 
labelled "fBm" (c.f. the supplementary material in [1^). 
However the solution of the equation for ftp is now known 
to be non-Gaussian [43[, and can be given in terms of 
Fox functions. Conversely the pdf of fBm is by definition 
pSl l46j Gaussian but with a variance which "stretches" 
with time as i^^. The correct kinetic equation for fBm 
must thus ^25] be local in time. It is shown in row 8, 
column 3 of Table 1. Given, to our knowledge, first by 
Wang & Lung [30|, it can be seen by trial solution to 
have a solution of the required form. 

The difference between ftp and fBm is striking, in that 
although both include temporal correlations, the kinetic 
equation for the ftp is non-local in time while that for 
fBm is local. This distinction disappears when we go 
to a Langevin description, where both processes explic- 
itly require fractional derivatives J23| • We are grateful to 
our referee for pointing out that the kinetic equation for 
fBm also corresponds to a transformation of time to t^^ 
in the ordinary diffusion equation for BWBm, which we 
may contrast with the fractional derivative in time in the 
kinetic equation for ftp. We remark that if one rescales 
BWBm with time, the resulting increments would not 
be stationary whereas fBm with the same kinetic equa- 
tion has stationary increments. This illustrates the point 
that fBm shares its kinetic equation with several other 
stochastic processes and so a full specification of the pro- 
cess thus requires more than the kinetic equation. 



D. Combining memory with infinite variance: 
"ambivalence" vs. Ifsm. 



1. "ambivalent processes" and the fully fractional 
continuous time random walk 



Questions similar to those in the previous section have 
been asked in the physics literature about the natural 
generalisation of the fractional time process to allow for 
both Levy distributions of jump lengths as well as power- 
law distributed waiting times. The resulting fully frac- 
tional kinetic equation is fractional in both in space and 



time: 



Direct differentiation of this pdf gives 



QaH ga 

-Q^Pap{x,t) = D Pap{x,t) 



(3) 



and was used by Brockmann at al. [12l | to exemplify the 
"ambivalent process" . Analogously with the ftp the solu- 
tion for this process is known [43| not to be a Levy-stable 
(or stretched Levy-stable) distribution but rather a con- 
volution of such distributions. 



2. A self-similar stable alternative to the a.p.: Linear 
fractional stable motion 

Analogously to the generalisation of the ftp to the am- 
bivalent process, there are several _ff-self similar Levy 
symetric a-stable processes, described in [43] ■ We here 
consider the simplest one, linear fractional stable motion 
(Ifsm), which generalises linear fractional Brownian mo- 
tion to the infinite variance case. We emphasise that Ifsm 
is, for example, not the fractional Levy motion referred 
to by \4§\. We can describe Ifsm through a stochastic 
integral: 



SniAaMM)^ KHAt-s)Z^ids) (4) 

J —CO 

where the memory kernel Kn.a is given by 

+ b.lit-sf-'/'^-i-s)"-'^'^] (5) 

Burnecki et al. [50| showed how mixed linear fractional 
stable motion can be obtained from Za (t) using the Lam- 
perti transformation [141, |5l|, a more general result which 
enables any self-similar process to be obtained from its 
corresponding stationary stochastic process. We are con- 
cerned here, however, simply with obtaining the kinetic 
equation. This can be found by direct differentiation of 
the characteristic function with respect to time (c.f. [52l|). 

As with the simpler stable processes the pdf pifsm of 
Ifsm can be expressed via the Fourier transform of the 
characteristic function (e.g. |13l.[53|) 



Pif. 



Akx 



exp(- 



-awe") 



(6) 



We see that the characteristic function: p = 
exp(— CT|fc|"i"'^) generalises the oLm case. Because a 
is no longer equal to 1/ H the effective width parameter 
now grows like t"^ . The characteristic function has the 
correct fBni limit, when a = 2, we see for fBm at any 
given t it is a Gaussian with width growing as t^^ . We 
can see that Ifsm is a general stable self-affine process by 
taking k' — kr^ which gives pifsm = t~^4'aix/t^), a 
stable distribution of index a and a prefactor ensuring 
_ff— selfsimilarity in time [28| . 



g-^Pifsm = -aHae"-' J e^'^^lkr exp(-a|A:rr^) 

(7) 
which, absorbing the constant a, and factors of 2tt into 
D can be recognised as 



^ Pifsm ^aHe^'-^D^pifs 



dt 



dx° 



(8) 



using the above definition of the Reisz derivative. Sur- 
prisingly the kinetic equation of Ifsm seems not have been 
given explicitly before in either the physics or mathemat- 
ics literature. Krishnamurthy et al. [28j quoted an equa- 
tion of motion for integrated activity in Ifsm. This has a 
more complicated structure presumably due to additional 
memory effects arising from the integration process. 

We note that aHt""-^ = dtt"'" ■ This factor arises 
because ([8]) could also be obtained from the space frac- 
tional diffusion equation ([21) by a simple transformation 
of the time variable: t is replaced by t"^. The appro- 
priate limits may be easily checked; in particular a = 2 
gives the kinetic equation of fBm. 

We also remark that Ifsm should be a special case of the 
nonlinear shot noise process studied by Eliazar & Klafter 
[54| which may allow further generalisation of the kinetic 
equation we have presented. 



III. LFSM AS A MODEL OF INTERMITTENT 
BURSTS 

Intermittency is a frequently observed property in 
complex systems, and can be studied within several 
paradigms. One such, of continuing interest, has been 
Bak et al's self-organized criticality (SOC), a key pos- 
tulate of which is that slowly driven, interaction domi- 
nated, thresholded dynamical systems will establish long 
ranged correlations via "avalanches" of spatiotemporal 
activity. The avalanches are found to obey power laws 
in size and duration. In consequence, many papers have 
sought to measure "bursts" of activity in natural time 
series. This has most typically been done by means of a 
fixed threshold. The duration t and size s of the bursts 
are then respectively defined as the interval between the 
i-th upward crossing time (ti) and the next downward 
crossing time (ti+i) of the threshold, and the integrated 
area above the threshold between these times. 

The search for SOC in the magnetosphere and solar 
wind has used this approach among others (e.g. [23.l34|). 
The similarity of observed burst size and duration distri- 
butions in solar wind and magnetospheric quantitities to 
those from models of turbulence and SOC led Freeman 
et al. [3J| and Watkins [23] to conjecture that, at least 
qualitatively, such behaviour might simply be an artifact 
of a self similar (or multifractal) time series, rather than 
unique to a given mechanism. In particular this was in 
distinction to the idea that one could use the presence 



the isoset thus scales as: 



* Simulation 
p = 2-H 



X 



0.1 0.2 0.3 0.4 



FIG. 1: Dependance of exponent /3 for pdf p(r) of a burst of 
duration r on H for simulated Ifsm in the fBm, a = 2 limit. 
The average of 7 trials was taken. 



where 



p{t) 



(3 = 2-H 



(9) 



(10) 



giving the same exponent 3/2 as for the first passage dis- 
tribution in the Brownian case. For symmetric processes 
this scaling is retained by the subset of the isoset that 
corresponds to burst "durations" . We expect this to be 
independent of the detailed nature of the model and so 
should, in particular, also apply to Ifsm. 

To establish the behaviour of burst "sizes" we note 
first that Kearney & Majumdar [37] considered the zero- 
drift Wiener Brownian motion (BWBm) case. Rather 
than their full analytic treatment, we recap their heuristic 
argument for a burst size (area) A defined using the first- 
passage time tf. This may then be adapted to burst sizes 
defined using isosets, and thence to Ifsm. They first noted 
that for BWBm the instantaneous value of the random 
walk y{t) ~ t^/^ for large t. Then, defining A by 



or absence of power laws in waiting times defined sim- 
ilarly to the above as a distinguishing feature between 
SOC and turbulence. One of the present authors has 
thus elsewhere ^^ advocated the testing of avalanche 
diagnostics using controllable self similar models. Simi- 
lar points have been made subsequently by Carbone et 
al. [33 for fBm and Bartolozzi ^] for the multifractal 
random walk. 



In this section we thus present a preliminary investi- 
gation of the ability of Ifsm to qualitatively mimic SOC 
signatures in data. As the kinetic equation we have de- 
rived is not unique to Ifsm and is insufficient to specify 
all its properties in what follows we have used a numer- 
ical simulation of the process S using the algorithm of 
Stoev & Taqqu [38] and analytic arguments, after those 
of Kearney & Majumdar [37|| to predict the scaling of 
the tail of the pdf of burst size s and duration t in 
Ifsm for large s,t. Rather than estimate the exponents 
from plots of the pdf or empirical complementary cdf 
we have elected to use maximum likelihood estimation, 
as imple mented in the algorithms of A. Clauset and co - 
workers ( |http://www.santafe.edu~aaronc/powerlaws7 1 • 
Detailed comparison with measured experimental expo- 
nents is not attempted at this stage and will be the sub- 
ject of future work. 



The expected behaviour of burst duration and sizes for 
Ifsm has, as far as we know, not been investigated. Deal- 
ing first with durations, we make use of the fact that for 
a fractal curve of self similarity exponent H and dimen- 
sion D = 2 — H the points {ti} have dimension 1 — H. 
In consequence the probability of crossings over a time 
interval r goes as t^~^ giving an inter-event probability 
scaling like t~^^~^K The pdf for inter-event intervals in 



y{t')dt' 



(11) 



,3/2 



the integration implies that large A scales as t / . Simple 
inversion of this expression implies that tf must scale as 
tf ^ A^l'^ . We independently have the standard result 



-3/2 



/ 



To 



for first passage time for BWBm: P{tf) ^ t 
get P{tf) as a function of A i.e. P{tf{Aj) one needs to 
insert the expression for tf as a function of A in above 
equation, and in addition will need a Jacobian. After 
these manipulations P{A) ~ ^""/^ f37J . 

In the zero-drift but non-Brownian case we will still 
argue that y{t) ~ t^ for large t. As our application uses 
the above mentioned isoset-based burst size s rather than 
those based on first passage times we define s by: 



t.+i 



y{t')dt' 



(12) 



The rest of the argument goes as before but using 
We find: 



where 



P{s) 



7 = -2/(l + ff) 



(13) 



(14) 



which we can check in the Brownian case where H — 1/2 
to retrieve P{s) ~ s~^/^. 

The same exponents, (3 and 7, but defining the bursts 
using a DFA-like moving average rather than a fixed 
threshold, were earlier investigated, for the fBm case 
only, by Carbone et al. [S^. The format of our figures for 
the fBm and Ifsm cases has been chosen to allow com- 
parison with theirs. They found numerically the same 
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FIG. 2: Dependance of exponent 7 for the pdf p(s) of a burst 
of size s on H for simulated Ifsm in fBm limit. Again an 
average of 7 trials was taken. 



FIG. 3: As Figure 1, dependence of duration exponent /3 on 
H, but 1 trial only. 



dependences of P and 7 on H, as we have in equations 
PH)) and (|14p above, which is intuitively reasonable with 
hindsight because the choice of fixed or running threshold 
should not change the asymptotic scaling behaviour. For 
a fixed threshold the burst size and duration exponents 
for fBm have also very recently been presented by [39|. 

Our numerics confirm that using the fixed threshold 
definition the expressions for f3 and 7 also describe fBm 
well, although the scatter, from a single trial in the case 
of each value of (3 and 7 shown in Figures 3 and 4, seems 
relatively high. We have reduced the scatter in figures 
1 and 2 by plotting the average of the exponents over 
a small number of trials (here 7) . The assumption that 
burst size s grows with duration r used in our heuristic 
derivation above can be seen to be reasonable for the 
fixed threshold, fBm case in Figure 5. 



Simulation 
-y.2/(1+H) 



0.1 0.2 



FIG. 4: As Figure 2, dependence of size exponent 7 on H, 
again 1 trial only. 



Perhaps more surprisingly the expressions also seem 
to hold reasonably well when the stability exponent is 
reduced, first to 1.8 (Figures 6 and 7) and then to 1.6 
(Figures 8 and 9). Again we note that these are averages 
of 7 trials in each case. By the a = 1 case presented in 
Figures 10 and 11, however, the expressions can be seen 
to fail. In this parameter regime, for any given H, they 
are seen to consistently underestimate both the burst ex- 
ponents. It has been suggested to us that this could be 
because y^ ceases to be a good estimate of character- 
istic displacement when the increments of the walk are 
very heavy tailed (S. Majumdar, Personal Communica- 
tion, 2006) but we have so far been unable to find a 
suitable alternative expression. 



IV. CONCLUSIONS 

In this paper we studied the question of whether one 
would expect the same equation to describe a time se- 
ries as an anomalous diffusive process. A codification of 
diffusion-like equations showed that a kinetic equation 
was "missing" from the literature; the one correspond- 
ing to Ifsm. We gave a simple derivation for it by direct 
differentiation of the well-known characteristic function 
of Ifsm. We then made a preliminary exploration of how 
Ifsm could model the "burst" sizes and durations previ- 
ously measured on magnetospheric and solar wind time 
series. We made simple scaling arguments building on 
a result of Kearney & Majumdar [371 to show how Ifsm 
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FIG. 5: Dependence of exponent ^p on Hin fBm case, tp 
captures growth of burst size s with duration r 



FIG. 7: As Figure 2 (burst size exponent 7 vs H, 7 trials), 
but for a = 1.8. 
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FIG. 6: As Figure 1 (burst duration exponent /3 vs H, 7 trials), 
but for a = 1.8. 



FIG. 8: As Figure 1 (burst duration exponent /3 vs H, 7 trials), 
but for a = 1.6. 



could be one candidate explanation for such "apparent 
SOC" behaviour and made preliminary comparison with 
numerics. These arguments fail when the tails of the pdf 
of increments become very heavy, and further work is 
needed on this topic. 

In future we also plan to consider other stochastic 
processes, both FARIMA (c.f. |2j|) and nonlinear shot 
noises, to allow generalisation of the above initial inves- 
tigations into burst size and duration. The prevalence of 
natural processes showing heavy tails and/or long ranged 
persistence suggests a relevance well beyond our initial 
area of application in space physics. 
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